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Simulating 
Using Neural 


summary 

models ° f Si T' at ' ng and synthesizi "8 substructures by computational neural network 

stress modelling Th/llT^ *** I” 8 f statlcally ‘"determinate beam, using both a 1-D and a 2-D plane 
g. he beam can be decomposed into two cantilevers with free-end loads By training 

beTlved I rf ma“ch“ U D Tet^ T^h ‘° diff<ire "‘ ' he ori S inal bea "> Pr°W<® can 

al^ritW match ; u P betw « en tw ° subsystems under compatible interface conditions. The genetic 

agreement Pr ° blem ' “ fo “" d to g “ d 


Introduction 

Many aerospace structures are composed of substructures that are geometrically simpler and 
u hose responses to loads are relatively easier to calculate. While a typical load-response relationship of 
the substructure can be established by utlizing neural network (NN) models trained from finite element 
na ysis data, how to mcorprate the NN models to get the characteristics of the whole system is a 
problem to be solved. We investigate the feasibility of synthesizing computational neural network models 
of aerospace substructures by studying some simple problems. For statically determinate structures in 
superposition can be applied and synthesizing of the NN models for the substructures 

The statically indeterminate structures, on the other hand, must be split into statically determi 
nate substructures first. When neural networks representing the substructures are trained, compatibility 
conditions needs to be enforced to obtain the unknown interface loads. But how to solve this match-up 

thiS b, PaPW a gl ° bal SearCh ,edmi<1Ue ' Le - the «“* X 

the interface quandties P “ “ ° pt ' m,zation P rocess of differences between 

\rvr ■ ^'” Ce ma j°r part of solving the statically determinate structures, i.e. specifying and trainning 

,1,3 also addressed in the statically indeterminate problems, only the latter will be addressed here. 

A 1-D Beam and Its Analytical Solutions 

solution^iffhpK^ indete . rminate bea ™ P roblem sh °™ ^ Fig. 1(a) has the following exact analytical 
B . ’ f h beam 18 uniform Wlth a bending rigidity of El. The interface moment and force at point 
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Substructure 2 


Fig. 1 A Statically Indeterminate Beam and its Decomposition 


M = Ifl aJ f,263 )( L l- L 2)-Pi a 2 (Lf-LiL 2 -2r,|)+P 2 6^f2 f,? + f„ I „_ / ,^ 

(Li+Li) 3 2L 

P _ ^i±L2}_(Pia 2 - P 2 b 2 )-2(P ]a 3 - P^b 3 ) 

(FTIir 


The displacements at point B: 


W B = f^ 3 ^(3^l+^)+P26 3 L?(L 1 +3 L 2 )-3P 1 a^ 1 Lj q 1 +L,l-3/>,fc^.?^fr.^r^ 

6£/(£ 1+ £, 2 )3 1 — ' a 

Or = Hi T2 fc3 - p l a3 )^i^+Pia 2 ^2(2Lf+L,L,- L^_p^2r |f _ r 2 +f|r . ,, r ^ 

6E/(^l+I 2 ) 3 ^ 


( 1 ) 


( 2 ) 


Simulation of the 1-D Beam by NN and GA 

Artificial Neural Networks (ANN), or simply Neural Networks (NN) are computational systems 

S rl' t ™ their data ■—* a " d restoring rnethod^ and learning 

S , T Y Slrable qUaIltiGS ° f NN ’ SUch 35 a universal approximator have enabled it beinf 
widely applied m engineering. For details of NN one can consult [1] S 

. m J°f mUl r! e ® tructare shown in Fi §- 1(a) by Neural Networks, the statically indeterminate 
are stat e n* substructures as shown in Fig. 1(b). These substructures (cantilevers here) 

F)are known mmate Pr ° Vlded the interface forces at P oint B (bending moment M and shear force 



Fig. 2 Cantilever Beam to be Simulated by a NN 


and two output variables: ’ g- • he network has three input variables 

{A/,F,a} =>NN=^ {w,d} (3) 

The neural network will be trained by using training data generated from the analytical solutions 


W B 


= -U/4-|F+ia 2 -ia 3 


N) 


Jb = M-\F-\a 2 
where all the quantities have been nondimesionalized: 

= o~ = EUb J7_ M - F_ a 

PL 3 ,e ' M ~PL’ F= p' a= L' 

1 hidden ^ayer'and lO^Md* ’T* *** & ML ^[ multi ' Iayer Perceptron, also called feed-forward N T N) with 
± uiauen layer ana iU hidden layer neuros. There arp ISO nf training . , . ' 1 

variations in M, F and a respectively S d ta SetS ’ a C0mbmatl0n of 6 x 6 x 5 

with L™nbtrg P S2dt bl tmW “ f ° Und ‘ ha ' the MATLAB Program for training feed-forward NX 
ing feed-forwari nZo* *1,^ c be 1 reSUltS ' 0ther programs ’ such “ (train- 

pe^formanc^wdi^t^^num^er^^input^vmdables^is 1 ^ ££ “* ^ 

force s ^‘™<=tures has been trained, the interface moment M and 

po,nt B can be obtained by enforcing the following compatibility and equilibrium conditions 

Mi = M 2 \8 BX = -8 B2 ;F l = —F 2 ;w B1 = w B2 . 

where the sub-index 1 and 2 indicate substructure 1 and 2 respectively. 

First the following approach was tried: 

{M,F} — ► {A/, F, a} =» NNl =» { w ,0} 

— > {w, 9, a} =» NN2 =* {M, F} 

^Tafa^ouLte o?l UFt* I' 11 ' Same Kheme 88 in Eq " (3 >. while NN2, having inputs of 

pair of vaLt (M F 1 ,M ’ “ ”7“ netW ° rk ‘° NN1 ' ThiS appr0ach did work since the 

Lial t FI ‘5’ f IJT not i C01 ; VCT 8 e - 11 was seen that even those iterations which started with an 
mtial {M, F} same as the exact solutions of Eqn. (4) and its reverse form 


M — 6 w B + 4 9 b — a 2 + a 3 
, F = 12 w B + 68 b — 3a 2 -I- 2a 3 


( 3 ) 


did n0t “ nv 3 because Eqn. (5) magnifies the error between the iterations enormously. 

iated a c -up problem in another way the following optimization process can be formu- 


where 


Miri\f pE(M, F) 

{M, F) - {]-[{w B2 - w m ) 2 + ( 0 B2 - 0 b ,) 2 ]}F 


( 6 ) 


( 7 ) 


• ' mis optimization prooiem can oe smvpn ^ 

stochastic methods based on the concepts of natural select, a, S omnms A-tAJ, tne aenvative-tree 

s “ *> “ - <» - ...““rxzr„z. -f-s 

used in the present study are described as follows: 2 ' Specificatl0ns of the GA 

related to T°r? 22 ' bit ^ S, ™S witl > <”» H-bit part 

or a gene. In each of the 11 bit string n, i tated t0 F 0ne blt of t,le string is called a chromosome 
1 b “ re ” the si «"' " ith » representing +, and 

non-zero values recognizable in this GA TZ mi mSgn “ Ude 10_< ' S ° ‘ he smallest 

strtng Sd1ntr F ^"S “* ^ * “ by a ^ 

higher the population. A member with a 

« /■ = wkKv where E is defined in Eqn “ ) ge " eratlon ' 11 b<> ^ 

vaiues. For the £35^1^^ haVe the beSt fitn " ss 

1W The^cro^oveMT^nK^hwhe^portioM o^the^ 1 " ° f V? ^ * * probability ° f 

In the first n' “* "-S. 1 -" ****■ °" * he P^formance of the GA. 

as functions of n, en were used and improvements were seen. ‘ Pm '"“ 


A 2-D Beam and Its Solution 

is shoJnln'KvTf 2°n t synthesizing substructures were used to solve a 2-D problem 
stresses “> ad P ' a “‘ ^ 

The shear and normal stress at A-A were assumed to have the following distribution 


where y = 


T xy — 


- a 2 {l - y 2 ),a x = b 0 + + b 2 y. 2 


( 3 ) 
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Fig. 3 A 2-D Beam under Parabolic Load 


Imagine that the beam is split at A-A and become 
Neural networks can be trained by taking FEM solutions 
combinations of a 2 ,6 0 ,6 x and b? as the inputs- 


two cantilevers with a shear and normal load, 
of the node displacements as the targets and 


{a 2 ,& 0 } =* NN 1 =» { u u vi at A-A, left cantilever}, 
{61,62} => NN 2 => { Ui , Vl at A-A, left cantilever}, 
{a 2 , 6 0 } => NN 3 => {u„ Vi at A-A, right cantilever}, 
{61, 6 2 } => NN 4 => {ui,Vi at A-A, right cantilever}. 
Then the interface stresses and deformation can be obtained by solving 


where 


6 0 , b u 62) 

E(a 2 , 6 0 , 61, 62) ~ ^ 2 l( u i C ~ wf c ) 2 -I- (v^ c - wf c ) 2 ]2 


A S a " d Tight ™ ntileveT ^tively. The GA used was 


Some Results 


For the 1 D problem the following cases were used: ^- = 0 8--*- = 0 8-& = n*-£z-n- 1 - 

Results are shown in Fig. 4 in which nart (»\ ;<? rt, Ll u u- ’ Ll r , Pl U ' 8, if ~ O.o ~ l.o. 
(T It — -i /u\ • , ,. ’ . n wnicn P art ( a ) 1S the search history of the GA for the first nnint 

X ^ ^ 1S ^ dlstributi °n of the population on the search plane when GA is terminated 

( ) e^comparison of simulaions and the exact solutions of M, and F and fdi is rtm . 1 

th * -‘“‘'I* * Micates that the nondtS^ ° f 
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(c) (d) 

Fig. 4 Calculation of forces and displacements at point B 

pared ^th FEM stow^nFigs TLt We^JaT T™^ 0 " “ interface A ' A »■»- 

the finite element analysis, and the mesh is shown in Fig 6 r «tangula.r element ([3]) to do 

Stress Distributions at A- A 
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Tauxy Simulated 
Sgmx Simulated 


Mesh before deformation 

Mesh after deformation by FEM (exaggerated) 

* * Deformed interface A-A by NN and GA simulation 
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Fig. 5 Stress Distributions at A-A 


Fig. 6 Mesh and Deformations 


6. Conclusion 

direct ZZSXSSZZ Can " 0t ^ ^ ^ SimP,£ “ ^ "* 

more “ “ U,d * ^ a PP« «o 
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